First-order super-radiant phase transitions in a 
multi-qubit-cavity system 

Chiu Fan Lee* and Neil F. Johnson^ 

Centre for Quantum Computation and Physics Department 
Clarendon Laboratory, Oxford University 
Parks Road, Oxford 0X1 3PU, U.K. 

February 1, 2008 



Abstract 

We predict the existence of novel first-order phase transitions in 
a general class of multi-qubit-cavity systems. Apart from atomic sys- 
tems, the associated super-radiant phase transition should be observ- 
able in a variety of solid-state experimental systems, including the 
technologically important case of interacting quantum dots coupled 
to an optical cavity mode. 
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Phase transitions in quantum systems are of great interest to the 
soHd-state and atomic communities 012], and have even caught the 
attention of the quantum information community in connection with 
entanglement [HI l^- Most of the focus within the sohd-state com- 
munity has been on phase transitions in electronic systems such as 
low-dimensional magnets [2] while in atomic physics there has been 
much interest in phase transitions in cold atom gases and in atoms 
coupled to a cavity. In particular, a second-order phase transition, 
from normal to superradiance, is known to arise in the Dicke model 
which considers N two-state atoms (i.e. 'spins' or 'qubits' cou- 
pled to an electromagnetic field (i.e. bosonic cavity mode) 
The Dicke model itself has been studied within the atomic physics 
community for fifty years, but has recently caught the attention of 
solid-state physicists working on arrays of quantum dots, Josephson 
junctions, and magnetoplasmas jH]. Its extension to quantum chaos 
0, quantum information ^Uj and other exactly solvable models has 
also been considered recently [TT] . 

Despite its obvious appeal, the Dicke model ignores interactions 
between the 'spins'. In atomic systems where each 'spin' is an atom, 
this is arguably an acceptable approximation if the atoms are neutral 
and the atom-atom separation d ^ a where a is the atomic diam- 
eter. However there are several reasons why this approximation is 
unlikely to be valid in typical solid-state systems. First, the 'spin' 
can be represented by any nanostructure (e.g. quantum dot) possess- 
ing two well-defined energy levels, yet such nanostructures are not 
typically neutral. Hence there will in general be a short-ranged (due 
to screening) electrostatic interaction between neighbouring nanos- 
tructures. Second, even if each nanostructure is neutral, the typical 
separation d between nanostructures in fabricated arrays is similar 
to the size of the individual nanostructures themselves (10^ — 10^ A). 
Hence neutral systems such as excitonic quantum dots will still have 
a significant interaction between nearest neighbors |12j . 

Motivated by the experimental relevance of 'spin-spin' interac- 
tions, we introduce and analyze a generalized Dicke Hamiltonian which 
is relevant to current experimental setups in both the solid-state and 
atomic communities We show that the presence of transverse 

spin-spin coupling terms, leads to novel first-order phase transitions 
associated with super-radiance in the bosonic cavity field. A techno- 
logically important example within the solid-state community would 
be an array of quantum dots coupled to an optical mode. This mode 
could arise from an optical cavity, or a defect mode in a photonic 
band gap material [E]- However we emphasize that the 'spins' 
may correspond to any two-level system, including superconducting 
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qubits and atoms [SJEj. The bosonic field is then any field to which 
the corresponding spins couple jHl EI- Apart from the experimen- 
tal prediction of novel phase transitions, our work also provides an 
interesting generalization of the well-known Dicke model. 

Just as the themodynamic results for the original Dicke model 
turned out to be valid for a wider class of Dicke-like Hamiltonians [7], 
the results we present are actually valid for a wider class of Hamilto- 
nian incorporating spin-spin and spin-boson interactions jl4j . How- 
ever for simplicity, we will focus here on a straightforward example: 

^ r A 1 

= + E |^(« + «^)^f + - j'^i ■ '^r+ij • (2) 

Following the discussions above, the experimental spin-spin interac- 
tions are likely to be short-ranged hence only nearest-neighbor in- 
teractions are included in H. The operators in Eqs. 1 and 2 have 
their usual, standard meanings. To solve for the thermodynamic 
properties of H, we introduce Glauber coherent states |a), which 
have the following properties [Tj: a\a) = a\a), {a\a^ = {a\a* , and 
/ ^^^(")^^™(") |a)(a| =1. In terms of this basis, the canonical parti- 
tion function can be written as: 

Z(iV,r) = E/ ^^'^"f"^^"^ sKa|e-^^|a)|s) (3) 

We adopt the following assumptions as in Ref. [7j: 

1. a/\/]V and a^/\/iV exist as ^ oo; 

2. limTv^oo limR^oo Z]f=o ^"^^^'^ can be interchanged 
We then find 

Z{N,T) = I ^e-'3|-l'(TVe-^^y (4) 
J vr 

where 

, ^ f ARe(a) ^ e z jY y \ 




= Eefe(«)(7i7.-^) (7) 
k=i ^ 
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with 



5(a) 



2jJl + ig{a))^+2g{a) 



'l^ARe(a) 



JVN 



(8) 
(9) 



Equation 6 follows from a rotation about the y-axis, while Eqs. [7|and 
IHlare derived in Ref. [2j for example. With H' diagonalized, Z{N,T) 
becomes 

[k=lnk=0 

= J ^g-/3(|ap-Efc«fc(R«(-))/2) I fj (l + e-^«^(^^("))) I (10) 
= _L^|d^e-/3(-^-E.«.(-)/2)|^(i + e-^«'=(-))| (11) 

where w = Re(a). Letting x = w/\/JV and writing J2k=i ^ lo^ dk, 
Z{N,T) becomes 




dx |e 



TV 



where 



and 



1 r 

= ^ dk |log 



cosh ( ^S,k{x] 



(12) 



+ log(2) (13) 



a(x) = 2 J Jl + {g{x)y + 25(x) cos k 



(14) 



Prom here on, we will omit the log(2) term in I{x) since it only con- 
tributes an overall factor to Z(N,T). 
Laplace's method now tells us that 



Z{N, T) oc max exp \N[-f5x^ + I{x)] ] 

— oo<x<oo ^. J 



(15) 



Denoting [— /3a;^+/(x)] by i7(x), we recall that the super-radiant phase 
corresponds to Q,{x) having its maximum at a non-zero x [7j- If there 
is no transverse field, i.e., if J = 0, and the temperature is fixed, 
then the maximum of Q.{x) will split continuously into two maxima 
symmetric about the origin as increases. Hence the process is a 
continuous phase transition. 
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However the case of non-zero J is qualitatively different from J = 
0. As a result of the frustration induced by the tranverse nearest- 
neighbor couplings, there are regions where the super-radiant phase 
transition becomes first-order. This phenomenon of first-order phase 
transitions is revealed by considering the functional shape of I{x). 
Numerical simulations show that I{x) can have one, two or three local 
maxima, as can be deduced from the behavior of in Fig. 1. 

Figure 2 shows the regions in which these three cases appear for 
fixed J and f3. These regions are determined by considering the sec- 
ond derivatives of I{x) for various A, J, e and j3. Figure 3 plots the 
maximizer of Vt{x) with A fixed at a value of 1.3. Referring to Fig. 
2, we see that in the region with three maxima, e.g. when e < 1, the 
non-zero local maximizers are dominant and hence the super-radiant 
state appears. As e increases, these two local maximizers converge to 
zero and the system is no longer super-radiant. This is no longer the 
case if J is increased slightly, e.g. to 0.56. In this case, Vt has a global 
maximum when e is small; however as e increases, the non-zero local 
maxima become dominant and as a result, a first-order phase transi- 
tion occurs. We note that the barriers between the wells are infinite 
in the thermodynamic limit, hence we expect that the sub-radiant 
state is metastable as e increases. This observation also suggests the 
phenomenon of hysteresis, which awaits experimental validation. 

The relevant parameter ranges J ~ 0.1 — 0.6, fi = 100, e ~ 1 
and A ~ 0.5 — 1.5 correspond to realistic scaled parameter values for 
both nanostructure-cavity and atom-cavity systems jl21 1131 ITC] . For 
example in semiconductor quantum dot systems, e is of order of the 
underlying bulk bandgap and hence can easily be chosen to range 
from 1.5eV in GaAs-like systems down to O.leV in narrow-gap semi- 
conductors. Meanwhile J is determined independently by the interdot 
separation which can be chosen to range from 10^ — 10^ Angstroms, 
hence J can be engineered to be of the order of ~ 10 — 10^ meV. The 
value of A can also be chosen independently, according to the strength 
of the coupling to the cavity mode. Both weak and strong coupling 
regimes are accessible using current nanotechnology. 

In conclusion, we have shown that the experimentally relevant 
spin-spin interaction transforms the nature of the super-radiant phase 
transition in the Dicke model. Our results highlight the importance 
of spin-spin coupling terms in spin-boson systems and open up many 
questions regarding the possible competition between the sub-radiant 
and super-radiant states in experimental systems. 

C.F.L. thanks University College for a Research Fellowship and 
NSERC (Canada) for partial financial support. N.F.J, thanks the 
LINK-DTI project, and is very grateful to Luis Quiroga and Alexandra 
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Figure 1: The function Q{x), which governs the phase transitions in the 
system, shown for various values of J and e with A fixed at 1.3. The y-axes 
are shifted for each curve so that all curves are zero at the origin. 
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Figure 2: Regimes of behavior of the function Q{x) which governs the phase 
transitions in the system. The red broken hues represent various values of 
J (J = 0.1,0.2,0.3,0.4,0.5,0.6 from bottom to top). Above the red broken 
hne and to the left of the blue solid line signifies the region where Q{x) has 
two local maxima. Below the red broken line and to the left of the blue solid 
line, signifies the region where Q{x) has three local maxima. Q{x) has only 
one maximum to the right of the blue solid line. P is taken to be 100. 
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Figure 3: Plot of the maximizer of the function Q{x), which governs the 
phase transitions in the system, as a function of J and e. A and (3 are fixed 
to 1.3 and 100 respectively. 
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